Correlation between muonic levels and nuclear structure in muonic atoms 
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A method that deals with the nucleons and the muon unitedly is employed to investigate the 
muonic lead, with which the correlation between the muon and nucleus can be studied distinctly. 
A "kink" appears in the muonic isotope shift at a neutron magic number where the nuclear shell 
structure plays a key role. This behavior may have very important implications for the experimen- 
tally probing the shell structure of the nuclei far away from the /3-stable line. We investigate the 
variations of the nuclear structure due to the interaction with the muon in the muonic atom and 
find that the nuclear structure remains basically unaltered. Therefore, the muon is a clean and 
reliable probe for studying the nuclear structure. In addition, a correction that the muon-induced 
slight change in the proton density distribution in turn shifts the muonic levels is investigated. This 
correction to muonic level is as important as the Lamb shift and high order vacuum polarization 
correction, but is larger than anomalous magnetic moment and electron shielding correction. 

PACS numbers: 36.10.Dr, 27.80.+W, 21.60. Jz,23.40.-s 

Keywords: muonic atoms, relativistie mean field, nuclear shell structure, isotope shift, muonic spectrum 



A muonic atom is an atom in which one of the electrons 
has been replaced by a negatively charged muon. Because 
of the very large mass of a muon compared with a electron 
and the correspondingly small Bohr radius, the muonic 
wave function has a large overlap with the nucleus. X-ray 
transition energies in muonic atoms are strongly affected 
by the size of the nuclei, and can be used efficiently to 
determine the nuclear charge distribution A de- 

tailed introduction about muonic atoms can be found in 
Ref. [S]. Moreover, the muonic atom tends to play an 
important role in investigation on other subjects. A re- 
cent study demonstrated that muonic atoms in strong 
laser fields can be used to dynamically gain structure in- 
formation on nuclear ground states [4;] . Masafumi Koike 
et al. proposed a new process of fi~e~ — > e~e~ in a 
muonic atom for a quest of charged lepton flavor viola- 
tion (This violation is known to be one of the important 
rare processes to search for new physics beyond the stan- 
dard model) in consideration of the fact that this process 
in a muonic atom has various significant advantages Q. 
An attractive means to improve the accuracy in the mea- 
surement of proton root-mean square radius is provided 
by muonic hydrogen Q . In addition, some potential syn- 
ergies of combining muons with radioactive nuclei may 
become a new tool to be used at future RIB facilities, as 
suggested in Ref. 0] ■ Therefore, investigation on muonic 
atom is meaningful and not limited to atomic physics. 

We extend the relativistie mean field (RMF) approach 
in this Letter to include the negatively charged muon, 
and then the correlation between the muon and nucleus, 
namely the effects of the nuclear structure on the muonic 
spectra as well as the influence of the muon on the nuclear 
structure, can be investigated distinctly. Within this cor- 



relation, we may find some new approach to probe the 
nuclear structure. Only one muon captured by the nu- 
cleus is discussed here. Because the motion of the bound 
muon is relativistie for high-Z atoms, it is necessary to 
treat it in relativistie framework. Nowadays the RMF 
theory has became a standard tool in low energy nuclear 
structure 8- 10] and the interacting Lagrangian density 
taking into account the muonic field is given by 

-^R.x ■ R"^ + \mlp, ■ p"" - gp'^-i^T ■ 

+i)^{i-l''dy - ■mp)iljf, - {-e)%lj^,Yi^^iK- (1) 

M, m„, and rup are the nucleon-, the cr-, the uj- 
and the p-meson masses, respectively. The muon mass of 
rUp = 105.6583668 MeV is taken from Ref. [ul- The 
nucleon field '0 interacts with the <J,uj,p meson fields 
CT, , pv and with the photon field A^, . The muonic field 
i>p interacts with the photon field A^, and A^, excludes the 
photon field produced by the muon itself. The field ten- 
sors for the vector meson are given as — d^uj\ — dxio^, 
and by similar expression for p meson and the photon. 
The self-coupling terms with coupling constants g2 and 
53 for the a meson are introduced which turned out to be 
crucial [13] ■ Since the muon couples only electromagnet- 
ically to nucleons, no additional parameters need to be 
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introduced, and no readjustment of the present param- 
eters is needed. Varying the effective Lagrangian, one 
can obtain the Dirac equation for nucleons and muons, 
and the Klein-Gordon equations for mesons. The calcu- 
lations are performed in coordinate space using a mesh 
size of 0.01 fm and different box sizes for different muonic 
states, and the paring correlations are accounted in the 
BCS formalism with an energy gap A obtained from the 
observed odd even mass differences for an open shell and 
A = for a closed shell. The NL3 and NLSH parameter 
sets are employed here. The NL3 parameter set has been 
used with enormous success in the description of a vari- 
ety of ground-state properties of spherical, deformed and 
exotic nuclei lla . Hj] , and NLSH is also a successful pa- 
rameter set 15]. The nucleons and muon are treated in 
the unified framework without any adjustable parameter 
and the finite size effect of the nucleus is automatically 
included so that both the muonic atom structure and 
nuclear structure can be investigated simultaneously. In 
contrast to earlier works about the muonic atoms based 
on Migdal theory [lB,[l3i the present method is relativis- 
tic. Also, our approach is more microscopic compared 
with the method that directly solves the Dirac equation 
using the two or three-parameter Fermi-type distribution 
of the nuclear charge [3 E^- In fact, it was a com- 
mon practice to fit the parameter sets for this nuclear 
charge distributions with the help of the experimental 
muonic levels [1^. In the RMF theory, the nucleons are 
treated as point particles, which is a drawback of RMF 
theory. As a matter of fact, the finite size effects of the 
nucleons cause some influence on the nuclear properties. 
The muonic spectrum is affected by the finite size ef- 
fects of the protons because the distributions of charges 
differ slightly from those of the protons. But in the self- 
consistent method, the proton density distributions are 
used as sources rather than the charge distributions. In 
our calculations, the Coulomb potential that the muon 
feels excludes the potential produced by the muon itself. 
In other words, the muon only interacts with the elec- 
trostatic potential generated by the protons. In fact, the 
mean field that the muon feels is a central potential field 
for a spherical nucleus, and the corresponding muonic 
energy i?^o is obtained by applying the RMF approach 
discussed above. 

Once the Dirac wave functions of the muon are ob- 
tained with the RMF method, the lowest order vacuum 
polarization correction Ae^, as shown in Fig. 1, can be 
calculated additionally. The details on this lowest order 
vacuum polarization correction can be found in Ref. [l9| . 
The ultimate calculated muonic levels = £"^0 + Ae^ 
in the muon-Pb systems are presented in Table I. As can 
be seen, the results of the RMF method taking into ac- 
count lowest order vacuum polarization corrections agree 
with the experimental data suggesting the effectiveness 
of the RMF to describe the muonic atoms. The lowest 
order vacuum polarization correction is relatively large 
for low energy levels, but the correction become smaller 
for high energy levels. The nuclear polarization that 
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FIG. 1: Lowest order vacuum polarization diagram. 



arises the muon-nucleus electrostatic interaction excites 
the nucleus into virtual excited states, is partly taken into 
account due to the self-consistent calculations, as will 
be shown below. Other corrections such as high order 
vacuum polarization correction. Lamb shift, anomalous 
magnetic moment, electron shielding and static hyperfine 
structure interaction, contribute little to the binding en- 
ergy of the muon atoms compared with the lowest order 
vacuum polarization correction [l^ and these effects are 
inessential for the following discussions, so they are not 
considered here. Note that the uncertainties of muonic 
levels especially lsi/2 levels in heavy nuclei are mainly 
ascribed to the nuclear structure, i.e., the proton den- 
sity distribution. The agreement between the calculated 
and experimental muonic levels indicates the reliability 
of the proton distribution from the RMF calculations. 
The lsi/2 levels with NL3 are in better agreement with 
the experimental data than those obtained with NLSH, 
possibly because that NL3 can provide a more accurate 
proton density distribution over the parameter set NLSH. 
The parametrization NL3 is more excellent because it 
cures some deficiencies of the NLSH to some extent, as 
suggested in Ref. [l^. In addition, the mean speed of 
the muon in lsi/2 orbit for the muon-^'^^Pb system is 
V = 0.31c according to our calculations, indicating the 
importance of relativistic effect. Thus the treating of the 
muon and nucleons in a unified relativistic framework is 
necessary. 

In Fig. 2(a), the nuclear isotope shifts {r'^)c — 
{r^)c(^°^Ph) for the Pb isotope chain are plotted as a 
function of neutron number, taking the nucleus ^'^^Pb as 
a reference. With the increase of the neutron number 
iV, the charge radius squared changes only slightly until 
it reaches the magic neutron number N = 126, but it 
increases more rapidly with the neutron number N goes 
beyond N = 126. In other words, a "kink" in the nu- 
clear isotope shift is found at iV = 126. We would like 
to mention that the non-relativistic calculations have not 
succeeded in offering this fact which originates from the 
nuclear shell effect. The muonic isotope shift (isotope 
shift of atomic spectrum) which is the shift of the muonic 



3 



TABLE I: The muonic spectrum obtained with the RMF method including the lowest order vacuum polarization correction for 
the muonic 204,206,208 pj^ rpj^g mean field part E^o and the lowest order vacuum polarization correction Ae^ are listed and the 
muonic levels are given by = E^o + Ae^. The experimental data 18] have been listed for comparison. All energies are in 
units of keV. 



A 


orbit 


^^^o(NL3) 


Ae^(NL3) 


£^o(NLSH) 


Ae^(NLSH) 


£m(NL3) 


i;^(NLSH) 


£^(Expt.) 


204 


lSl/2 


-10549.41 


-66.35 


-10561.24 


-66.47 


-10615.76 


-10627.72 


-10614.88 ± 0.41 


204 


2pi/2 


-4781.17 


-31.26 


-4783.41 


-31.30 


-4812.43 


-4814.71 


-4818.62 ± 0.08 


204 


2p3/2 


-4597.97 


-28.70 


-4599.67 


-28.73 


-4626.67 


-4628.40 


-4632.72 ± 0.06 


204 


2Si/2 


-3583.87 


-18.26 


-3585.90 


-18.28 


-3602.13 


-3604.18 


-3604.14 ± 0.26 


204 


3(^3/2 


-2161.81 


-9.40 


-2161.86 


-9.40 


-2171.21 


-2171.26 


-2173.03 ±0.05 


204 


3p3/2 


-2082.71 


-9.20 


-2083.28 


-9.21 


-2091.90 


-2092.49 


-2092.41 ± 0.33 


206 


lSl/2 


-10540.54 


-66.25 


-10551.51 


-66.37 


-10606.80 


-10617.87 


-10604.34 ± 0.41 


206 


2pi/2 


-4779.84 


-31.23 


-4781.97 


-31.27 


-4811.07 


-4813.24 


-4817.19 ±0.06 


206 


2p3/2 


-4597.03 


-28.68 


-4598.65 


-28.71 


-4625.71 


-4627.37 


-4631.50 ±0.05 


206 


2Si/2 


-3582.23 


-18.24 


-3584.08 


-18.26 


-3600.47 


-3602.34 


-3601.54 ± 0.34 


206 


3(^3/2 


-2161.78 


-9.40 


-2161.84 


-9.40 


-2171.18 


-2171.23 


-2173.02 ±0.05 


206 


3P3/2 


-2082.38 


-9.19 


-2082.93 


-9.20 


-2091.57 


-2092.13 


-2092.01 ± 0.39 


208 


lSl/2 


-10532.13 


-66.16 


-10542.23 


-66.26 


-10598.29 


-10608.49 


-10593.13 ±0.41 


208 


2pi/2 


-4778.56 


-31.21 


-4780.59 


-31.25 


-4809.77 


-4811.83 


-4815.14 ± 0.06 


208 


2p3/2 


-4596.13 


-28.66 


-4597.68 


-28.69 


-4624.79 


-4626.37 


-4630.25 ± 0.05 


208 


2Si/2 


-3580.66 


-18.22 


-3582.35 


-18.24 


-3598.88 


-3600.58 


-3599.75 ±0.17 


208 


3(^3/2 


-2161.76 


-9.40 


-2161.81 


-9.40 


-2171.15 


-2171.21 


-2173.01 ±0.05 


208 


3p3/2 


-2082.07 


-9.19 


-2082.59 


-9.20 


-2091.25 


-2091.79 


-2092.20 ± 0.32 
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FIG. 2: (a) Nuclear isotope shifts of the Pb isotope chain and 
(b) muonic isotope shifts of lsi/2 state in muonic Pb. The 
^'^*Pb is taken as a reference nucleus and the experimental 
data are taken from Ref. [lal. 



energy with increasing neutron number, includes the field 
shift coming from the change of spacial distribution of the 
charges in the nucleus and mass shift originating from the 
change of the mass of the nucleus. For a heavy nucleus, 
the field shift is major and much larger than the mass 
shift. In heavy elements such as lead, the mass effect 
is negligible and the field effect roughly accounts for the 
observed shifts. Fig. 2(b) displays the muonic isotope 
shift taking the muonic energy in ^"^Pb as a reference. 
Hence computational errors along with various correc- 
tions cancel to a high degree. This muonic isotope shift 
can be attributed to the nuclear isotope shift. The larger 
the nuclear charge radius is, the wider the charge dis- 
tributes since the total charge is invariant, and hence a 
high muonic energy. A "kink" in the muonic isotope shift 
is also found at = 126, which stems from the neu- 
tron shell effect in the final analysis. For other isotope 
chains, the "kink" can be also found at magic number, 
and we do not discuss them in detail here. This is cer- 
tainly a good news for the nuclear structure study. The 
"kink" in the muonic isotope shifts may be taken as a use- 
ful signal to experimentally probe the shell structure of 
the nuclei far away from the /3-stable line for which novel 
nuclear structure effects may exist. This is one of the 
main conclusions we draw in this work. Strasser et al. 
has proposed the cold hydrogen film method to extend 
muonic atom spectroscopy to the use of nuclear beams 
including radioactive isotope beams to produce radioac- 
tive muonic atom in the future (20| . This would allow 
studies of unstable nuclei by means of the muonic X-ray 
method at facilities where both /z~ and radioactive iso- 
tope beams would be available. The electron scattering 
is not going to be easy to be applied for exotic nuclei 
with very low beam intensities when one measures the 
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scattering cross sections to obtain the knowledge about 
the nuclear structure. But for muonic atoms, one mea- 
sures the X-ray energies, which is much easier and X-ray 
energy can be measured with a high accuracy. For the 
other approach to probe the nuclear shell structure, the 
first excited 2+ state, can be carried out using Coulomb 
excitation J2l\ for these exotic nuclei. Therefore, as two 
independent methods, the approaches of Coulomb exci- 
tation and muonic atom spectroscopy can complement 
each other. 



TABLE II: The correction for muonic Isip orbit due to the 
muon-induced slight change of the proton density distribu- 
tions in muonic Pb. 



Nucleus 


NL3 


NLSH 




-1.5 keV 


-1.2 keV 


206pb 


-1.5 keV 


-1.2 keV 




-1.5 keV 


-1.2 keV 


210pb 


-1.5 keV 


-1.2 keV 



Apart from obtaining information on the muonic spec- 
trum, the influence of the muon on nuclear structure can 
also be investigated. It is found that the nuclear struc- 
ture of the Pb isotope remains basically unchanged in the 
presence of the muon in lsi/2 orbit. The single particle 
level spacing of neutrons and protons is altered by only 
several or some dozen keV. The nucleon densities are re- 
duced by about 0.1% at the edge of the nucleus (~7 fm) 
but enhanced in the interior of nucleus by about 0.1%. In 
other words, the protons and neutrons move quite slightly 
towards the core as a result of the bound muon, and the 
nuclear rms radii are reduced only by about 0.02%. For 
the calculation of these changes in value, the systematic 
computational errors here could be canceled to a large 
extent, so that the results are reliable. The bound muon 
in other orbits affect the nuclear structure much more 
weakly than that in the lsi/2 orbit. The muon thus can 
be taken as a clean and reliable probe to extract infor- 
mation on nuclear structure since this probe can hardly 
change the nuclear structure. It is reliable to yield in- 
formation on nuclear charge distributions and rms radii 
through muonic X-ray in experiments. That is to say, 
the previous experiments on the nuclear charge distri- 
butions and charge rms radii within the muonic atoms 
were theoretically confirmed to be reasonable and reli- 
able. Moreover, the fact that the nuclear structure is 
only slightly changed by the muon iniplies that the lin- 
ear response theory used in Ref. [H, [13] is an excellent 
approximation. Although the change of the proton den- 
sity distribution due to the bound muon is quite slight, it 
in turn lowers the muonic levels especially for lsi/2 orbit, 
which is a part of nuclear polarization since only virtual 
excitation of = 0+ states is taken into account here. 
This correction for muonic Pb is presented in Table II, 
which is more important than the anomalous magnetic 
moment and electron shielding, and can be compared to 



Lamb shift and high order vacuum polarization correc- 
tions (For the lsi/2 orbit in muonic ^''^Pb, the correc- 
tions to muonic energy levels from anomalous magnetic 
moment. Lamb shift and high order vacuum polarization 
are 0.445 keV ,2g], 2.302 keV [H and -1.75 keV 
respectively. The electron shielding correction is even 
less, for example, -4.6 eV for the lsi/2 orbit in muonic 
2°^Bi [M]). However, this effect tended to be neglected in 
many previous investigations which regarded the charge 
distribution as an invariant when a muon is captured. 
Yet it is not easy to be studied with the usual method 
but can be automatically included in our approach since 
the muon and nucleons are treated in a unified theoret- 
ical framework. Attention should be paid attention to 
this correction since the experimental measurements on 
muonic levels can achieve this precision. 

In summary, the muonic atoms have been investigated 
in the framework of many-body theory. The main con- 
clusions are summarized as follows: (1) The relativistic 
mean field theory for nuclei has been extended to inves- 
tigate the muonic atoms of lead isotopes by taking the 
nucleons and a bound muon as a system and the theo- 
retical muonic levels agree with the experimental data. 
Therefore, a direct link has been established between the 
nuclear physics and atomic physics. (2) The muonic iso- 
topes shift stems from nuclear isotope shift for heavy 
nuclei. Just like the "kink" in nuclear isotope shift, a 
"kink" also appears in the muonic isotope shift at the 
neutron magic number. Perhaps it gives us a new ap- 
proach to experimentally probe the shell structure of the 
nuclei far from the /3-stable line via this "kink" in muonic 
isotope shift. (3) The muon-induced changes of nuclear 
structure are quite slight, which indicates the muon can 
be taken as a clean and reliable probe to extract infor- 
mation on nuclear structure, and theoretically confirms 
that the previous experiments on the nuclear charge dis- 
tributions and charge rms radii within the muonic atoms 
were reasonable and reliable. (4) As a part of nuclear 
polarization, the correction that the muon-induced slight 
change in the proton density distribution in turn shifts 
the muonic levels is studied, which should be paid atten- 
tion to because the contribution of it to muonic levels is 
more important than anomalous magnetic moment and 
electron shielding, and can be comparable to the Lamb 
shift and high order vacuum polarization correction. 
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